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Abstract. Consider the set of equations describing Oldroyd-B fluids in an exterior domain. It is 
shown that this set of equations admits a unique, global solution in a certain function space provided 
the initial data, but not necessarily the coupling constant, is small enough. 



1. Introduction 

The theory of Oldroyd-B fluids recently gained quite some attention. This type of fluids is described by 
the following set of equations 

Re(d t u + (u ■ V)w) - (1 - a)Au + Vp = divr inftx(0,T), 

divu = in ft x (0,T), 

We(d t T + (u- V)T + 5 a (T,Vu))+T = 2aD(u) inftx(0,T), 

u = on 5ft x (0,T), 

u(0) = uq in ft, 

t(0) = tq in ft, 



(1.1) 



where ft C K 3 is a domain, T £ (0, oo] and Re and We denote the Reynolds and Weissenberg num- 
ber of the fluid, respectively. Moreover, the term g a is given by j„(t,Vu) := tW(u) — W(u)t — 
a (D(u)t + tD(u)) for some a £ [—1,1], where D(u) = l/2(Vw + (Vu) T ) denotes the deformation 
tensor and W(u) = l/2(Vit — (Vit) T ) the vorticity tensor, respectively. 

This set of equations originally was introduced by J.G. Oldroyd |17j . whose intention it was to describe 
mathematically viscoelastic effects of certain types of fluids. 

The study of the above set of equations started by a pioneering paper by Guillope and Saut in 1990, 
see [5], who proved for the situation of bounded domains ft C R 3 , the existence of a local, strong solution 
to equation (jl.ip in suitable Sobolev spaces iJ s (ft). Moreover, this solution exists on [0, oo) provided the 
data as well as the coupling constant a between the two equations are sufficiently small. For extensions 
to this results to the L p -setting, see the work of Fernandez-Cara, Guillen and Ortega [2]. 

The existence of global weak solutions in the case of ft = R™ was proved by Lions and Masmoudi in 
[15] for a — 0. For extensions of this result to scaling invariant spaces of the form L^ c ([0, T); H S (M. 3 )) for 
s > 3/2, we refer to the work of Chemin and Masmoudi [T]. An improvement of the Chemin-Masmoudi 
blow-up criterion was presented recently by Lei, Masmoudi and Zhou in |11) . 

The situation of infinite Weissenberg numbers, mainly using the Lagrangian setting, was considered 
for ft = R 3 or for bounded domains ft C M. 3 by Lin, Liu and Zhang in [13], Lei, Liu and Zhou in [TO] 
and Lin and Zhang in |14| . Further results, describing in particular the two dimensional situation, can 
be found in in [12"]. [S] and 0. 
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The situation of exterior domains was considered first in [6]. There, the existence of a unique, global 
solution defined in certain function spaces was proved provided the initial data as well as the coupling 
constant a are small enough. 

In this paper we improve the main result given in [5] to the situation of non-small coupling coefficients 
a. Indeed, given an exterior domain fl C R 3 with smooth boundary, we prove the existence of a unique, 
global solution to equation (jl.ljl only under the assumption that the intial data uq and tq are sufficiently 
small in their natural norms. The idea of our approach is not completely new and goes back to the work 
of Molinet and Talhouk [IB]. The main idea is to take the divergence in the third equation in (| 1 . 1 1) and 
to control Pdivr in the i/^-norm by the corresponding term in L 2 and by curldivr, also measured in 
the L 2 -norm. For doing this, we need in particular to extend the well known estimate for the i? -norm 
of a function u by the L 2 -norms of u, div it, curlu and by u ■ v in the iJ^-norm for bounded domains 
to the situation of unbounded domains. Note that the main difficulty in the case of exterior domains 
is due to the failure of Poincare's inequality in this situation. Hence, we must treat the lower order 
terms of u and dtu in a new way and some new higher order energy estimates for u and d t u need to be 
developed as well. 

These higher order energy estimates then imply that the local solution is satisfying a certain differ- 
ential inequality which in turn implies that the local solution of (jl.l|l can be extended for all positive 
times. 

It is worthwhile to mention that also our approach to obtain local solutions to {HI)) seems to be new 
and quite different from the ones known in the literature. In fact, we base our local existence argument 
not on Schauder's fixed point theorem as it has been done in almost all of the existing works, but on a 
variant of Banach's contraction principle, avoiding hereby unnecessary uniqueness arguments. 

Some words about the derivation of equation (|I are in order. In fact, incompressible fluids are 
subject to the following system of equations 



(1.2) 



dtu+(u-V)u = diver, 
divu = 0, 



where u, a are the velocity and stress tensor, respectively. Moreover, a can be decomposed into —plal+r, 
where p denotes the pressure of the fluid and r the tangential part of the stress tensor, respectively. For 
the Oldroyd-B model, r is given by the relation 

(1.3) T + X 1 ^ = 2Jd(u) + X 2 D « D(u) 



Dt \ Dt 

where ^^ denotes the "objective derivative" 

^ =d t r+(u- V)r + 3a (r,V W ), 

and g a (r,Vw) := tW(u)—W(u)t— a (D(u)t + tD(u)) for some a 6 [—1,1]. Here D(u) is the deformation 
tensor defined as above and W{u) denotes the vorticity tensor. The parameters Ai > A2 > denote the 
relaxation and retardation time, respectively. 

The tangential part of the stress tensor r can be decomposed as r = rjv + r e where = 2n^D(u) 
corresponds to the Newtonian part and r e to the purely elastic part. Here 77 denotes the fluid viscosity. 
Denoting r e = r, the above equations can be rewritten as 



(1.4) 



dtu + (it • V)it — 77(1 — a)Au + Vp = divr, 

divu = 0, 



with the retardation parameter a := 1 — ^ G (0, 1). 

Using dimensionless variables, Oldroyd-B fluids may be thus described by the equations f| 1 . 1 1) . 
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In order to formulate our main result, let A be the Stokes operator in L 2 (£l) denned by 

Au := -PAu for all u G D(A) := H 2 (n) n #o(fi) n iv(fi), 

where the space L^(fl) is defined precisely later in Section 2. Moreover, we set V := i?o(0) n £^(f2). 
Our main results then reads as follows. 

Theorem 1.1. Let fl C R 3 6e an exterior domain with boundary dfl of class C 3 and let J = [0,oo). 
Then there exists eo > such that if 

(1-5) I^oIId(A) + ||7"o||ij3(n) < £o, 

then equation hl.l}) admits a unique, global strong solution (u,p,r) for all t G J satisfying 

u G C 6 (J; £>(A)) with Vu G L 2 ( J; # 2 ) and v! G L 2 (J; V) n C 6 (J; L 2 (ft)), 
Vp G L 2 (J;H 1 (fl))nL 00 (J;H 1 (Q,)), 
r G C b (J;H 2 (n))r\L 2 (J;H 2 (n)) with t' E C b (J; H 1 ^)) D L 2 (J; L 2 (n)). 

2. Preliminaries 

We start this section by recalling a higher order elliptic regularity estimate for the stationary Stokes 
equations. A proof can be found for example in [3], Theorem V 4.7. In the following, we denote by 
H k (il) the homogeneous Sobolev spaces of order k. 

Lemma 2.1. Let m G {0, 1} and n > 3. Assume that fi C K™ is an exterior domain with boundary of 
class C m+2 and g G H m (Q). Then the equation 

—Au + Vp = g in Q, 
divu = in 17, 
u = on dfl 

admits a solution (u,p) G H m+2 (Q) x H m+1 (Q,) which is unique provided Vit G L 2 (Q). In this case, 
there exists a constants C > such that 

(2.1) \\V 2 u\\ H ™ + \\Vp\\ H ™ < C(\\g\\ H ™ + \\Vu\\ L 2). 

The following variant of Lemma 12.11 will be important in the sequel. We first introduce the spaces 
G 2 (Q) := {u G L 2 (fl) : u = Vtt for some tt G #^(0)}, 

and 

L 2 (ft) := {u G C* c °°(^) : div u = in fi} - l|Ml ^ 
Then L 2 (J7) can be decomposed into 

L 2 (ft) = L^ft) © G 2 (0), 

and there exists a unique projection P : L 2 (Vl) — > L 2 (ft) with G2(f2) as its null space. P is called the 
Helmholtz projection. It is well known that 

L 2 (ft) = {ue L 2 {n) :dWu = 0,u-u\ 9U = 0}, 

where v is the exterior normal to dfl. 

Corollary 2.2. Assume that the assumptions of Lemma \2.1\ hold. Then there exists a constant C > 
such that the solution (u,p) of equation (|2.1[) satisfies 

(2.2) \\V 2 u\\ Hm <C(||Pg|| Hm + ||Vw|| La )- 
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Proof. Applying the Helmholtz decomposition to g and noting that P acts as a bounded operator on 
if" 1 (f2) yields g = Pg + for some V</s <E H m (ft). Hence, the first line of the above Stokes equation 
may be rewritten as — Ait + Vj> = Pg, where p = p — ip. Thus, estimate (|2.2[) follows immediately from 
eestimate (|2.1|) . □ 

Remark 2.3. Rewriting the above Stokes equation (|2.1|) as Au = Pg, we may replace the term Pg in 
(|2~2| by Au. 



We next recall a well known estimates for the if^-norm of a function defined on a bounded domain 
CI c M 3 . 

Lemma 2.4 f[3]). Lei f2 C M 3 be a bounded domain with boundary dCl of class C 2 . Then there exists 
a constant C > such that 

IM|iZ»»(n) < C (jMU^Q) + \\divu\\ L 2 {n) + ||cur/u|| L 2 (n ) + ||u ■ Hl ff 4( 9 n)) j u € if 1 ^) 3 . 



The following variant of Lemma |2 .41 concerns exterior domains. More precisely, the following propo- 
sition holds true. 

Proposition 2.5. Let fl C M 3 6e an exterior domain with boundary dCl of class C 2 . Then there exists 
a constant C > such that 

IMIifi(n) < C (JMIl 2 (q) + ||d«iHU 2 (Q) + ||cwrZu||L2(n) + ll u 1 Hl^i^ 

), ^^(O) 3 . 

Proof. We establish first a similar inequality for functions defined on K 3 . Assuming, for the time being, 
that u is smooth, we obtain for the j-th component u J of u 



Hence, 



and 



Au 3 ' = 9i(diU^ ' -d j u i )+d j divu = ^ <%(curlw) lJ ' + djdivi 

l<i<3 l<i<3 

u j = - A) _1 (curlu) ij ' - e?j(— A) -1 divu, 



<9 fe u j = - Sfc9i(-A) _x (curlu) <:7 ' - 9,-9k(— A) -1 divu 

l<i<3 

= - (cur ht) tJ ' — RjRkdivu, 

l<i<3 

where i£j := dj(— A)~i denotes the j-th Riesz transforms. Classical results on the boundedness of the 
Riesz transforms imply 

(2.3) llVulUa < C(||divu|| L 2 + ||curl«|| £ ») 

for some C > 0. By density, this estimate transfers to functions u belonging to H 1 (M. 3 ). 

Next, our aim is to combine the above estimate with the corresponding one in bounded domains in 
order to obtain the assertion for exterior domains. To this end, let R > such that fl c C Br(0) := {x £ 
R 3 : \x\ < R} and set 

D := n n B R+S (Q), Ti := and T 2 ■= dD\T±. 
Next, we choose a cut-off function <fr £ (B_r+3(0)) such that < <fr < 1 and 

[0, \x\ >R + 2. 

and decompose u as 

u = (f>u + (1 — <p)u =: + ti 2 . 
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Clearly, u\ € H 1 (D) and u\ = near IV Moreover, U2 € H^fl) and after zero extension (still denoted 
by U2), we can regard U2 as an element in i/ 1 (R 3 ). Now it follows from ()2.4j) and (|2.3p that 

|| Vzi||i3(o) < ||Vui|| L 2 (£)) + ||Vu 2 ||l2(r3) 

< C[\\<j>u\\ L 2 (D) + \\div(^m)\\ L2(D) + ||curl(0u)|| L 2 (D) + \\((jm) ■ v\\ H i 

+ 11 (<M ■ HI H £ (r2) + lldiv ((1 - 0)w)|| L 2 (R3) + ||curl((l - 0)u)|| i2(R3) ] 

< ' ^il ff i (ri) + \\4>u\\ L 2 (D) + ||V</> • U||z,a(£>) + ||V</>(g> U- It (8 V<f>\\ L 2( D ) 

+ ||^divu||i2(u) + ||(1 - 0)divu|| L 2 (R 3 ) + \\<pcm\u\\ L 2 (D) + ||(1 - ^curluU^^)] 

< ' + IMU^n) + ||divu|| L 2 (n) + ||curlu|| L 2( n )], 

and the estimate ()2.5p follows immediately. □ 

In order to construct a local solution to we study first two linearized equations; the first one 

for the velocity u and the second one for the tangential part of the stress tensor r, respectively. 
First, given T > 0, we recall some results on the Stokes equation 

' d t u - Au + Vp = / in Q x (0, T), 
divu = in Q, x (0,T), 

u = on an x (0,T), 

u(0) = tio in fl, 
where / is a given external force. 

Proposition 2.6 f|18|). Let fl C K 3 be an exterior domain with boundary dfl of class C 3 . Assume 
that f G L^TjIf 1 ^)),/' G ^(Q.T-.H- 1 ^)) and u Q G ff 2 (0) n V(fi). Tften i/iere exists a imzgue 
solution (li,p) of equation \2.1$ satisfying 

u G L 2 (0,r;F 3 (n))nC([0,r];iJ 2 (f7)ny(^)),u' G L 2 (0,T;F(O))nC*([0,T];L 2 (O)), 
p G L 2 (0,T;Hl c (n)). 
Moreover, there exists a constant C > smc/i t/iai 

ll«llia(fl»)ni«(fl»nv) + \W\\h(v)nL°°(Li) + IIVplH^ < C[\\u \\ 2 H 2 + \\f(0)\\v + \\f\\ L ^) 



(2.4) 



Next, consider the transport equation 

J We (dtT + (v ■ V)t) + t = 2<xD(v) - We ff (r,Vu), infix(0,T), 

\ r(0) = T , in 0, 

where t> is a given velocity field. 



(2.5) 



Proposition 2.7 ([18]). Lei SI C R 3 6e an exterior domain with boundary dfl of class C 3 . Assume that 
v G L 1 (0, T; _ff 3 (il) fl V(f2)) and tq G if 2 (f2). XTien t/iere exists a unique solution of equation \2. 5]) and 
a constant C > such that 

2a 

\\t\\l°°(p,T;H*((i)) < (\\to\\h* + ■£^r-)exp(C\\v\\ L i iH i } ). 

If in addition, v G C([0, T];H 2 (Q) n 7(0)), tfien r' G C([0,T\; H 1 ^)) and 

1 2a 

lk'llL~(o,T;i/i) < C(\\v\\ L ~ {H 2) + ^r^)(ll r o!lff 2 + ^T^-)exp(C||u|| L i ( H3)). 

The assertions of Propositions 12.61 and 12.71 are stated e.g. in [TH| even for a more general class of 
domains, however, without giving a proof. 
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3. Existence and Uniqueness of a Local Solution 

In this section we prove that the system (jl.lj) possesses a unique, local solution provided the intial data 
are smooth enough. More precisely, the following result holds true. 

Proposition 3.1. Assume that Q C R 3 is an exterior domain with boundary dil of class C 3 . Let 
uq G D{A) and To G H 2 (fl). Then there exist > and functions 

u G L 2 (0,T*;# 3 (fi))nC([0,T*];D(A)) with u' G L 2 (0, T*, V) H C([0, T*]; L 2 (0)), 
p 6 L 2 (0,T;H? oc (n)) with Vp G L 2 (0, T*; ff 1 ^)), 
r G C([0,T i ];if 2 (fi)) wf/i r' G C([0,T Jt ],ff 1 (n)) 
s«c/j i/iai (u,p, r) is £/ie unique solution to equation on (0,7^). 

Let us begin the proof of Proposition 13.11 with the following variant of Banach's fixed point theorem. 
For a proof, we refer e.g. to [7]. 

Lemma 3.2 ([7j). Let X be a reflexive Banach space or let X have a separable pre-dual. Let K be a 
convex, closed and bounded subset of X and assume that X is embedded into a Banach space Y . Let 
$ : X — > X map K into K and assume there exists q < 1 such that 

\\$(x) - <f>(x)\\ Y < q\\x-y\W, x,y G K. 
Then there exists a unique fixed point of <5> in K. 

Our proof of Proposition 13.11 relies on a combination of Propositions 12.61 and 12.71 with Lemma 13.21 To 
this end, consider for T > the following function spaces 

Ei := L 2 (0,T;H 3 (n))nL°°(0,T;H 2 (n)nV), 
E 2 := L 2 (0,T;V(n))nL°°(Q,T;Ll(n)), 
Fx := L°°(0,T;H 2 (n)), 
F 2 := L°°(0,T;ir 1 (fi))), 
and for B 1: B 2 > define the set K(T) by 

K{T) := {(v,9) e£ix Fx,v' G E 2 ,6' G F 2 ,v(0) = u ,9(0) = r and 
\\v\\ 2 El + \\v'\\ 2 E2 < B X , \\9\\ Fl < B u \\0'\\ F2 < B 2 }. 
Next, given (v, 9) G K(T), we define the mapping 

$M) := (u,t), 

where (u, r) is defined to be the unique solution of the corresponding linearized problem of II) 

Re d t u + (1 - a)Au = -Pdiv (v <E> v) + Pdiv 9 in X (0, T), 

We (StT + (u • V)r) + r = 2aD(u) - We g a (r, V«) inOx(0,T), 

(3.1) <( u = on dfl x (0,T), 

m(0) = wo in O, 

r(0) = To in 0, 

where A denotes the Stokes operator defined as in Section 1. It follows from Proposition 12.61 and 12.71 
that for appropriate choices of B\ and B 2 , there exists T\ > such that §(K(Ti)) C K(Tx). 

Next, we will prove that there exists T* G (0, Ti] such that $ is contractive on Y(T+), where Y(T) is 
defined by 

Y(T) := {(v,9) G L°°(0,T;L 2 (f2)) x L°°(0, T; L 2 (Q)), Vv G L 2 (0, T; L 2 (Q))} . 
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Indeed, for (vi,6i) G K(T\) let (ui,n) = §>(vi,6i) for i = 1,2. Moreover, we set u = u\ — u 2 and 
f = ti — 72 . Then (u, f ) satisfies the equation 



(3.2) 



Re d t u + (1 - a)Au = -Pdiv (v ® + v 2 ® u) + Pdiv6> in ft x (0,T), 
We Stf + f = 2aD(u) - We [(u • V)n + (u 2 ■ V)r+ 

3„(f,Vwi)+fl (r 2! V«)] inOx(0,T), 

u = on 9f2 X (0,T), 

u(0) = in 0, 

f(0) = in ft. 



Taking the L 2 inner product of (|3.2p i with u, we obtain 
~(Re \\u\\ 2 L2 ) + (1 - a)\\Vu\\ 2 L2 = Re(v®vi+v 2 ®v\Vu)-(9\Vu) 

< ^^IIVull^ + ^-(lkll£~ + IMIi~)Nl£* 

2 1 — a 



l-a 



L 2 



1 — OL ,,„ , 

< -5-11 V«|| 



CRe 2 
1 -a 



Mllr* + H«2ll!rOHi» 



1 — a 



11*1 



Consequently, 

(3.3) |(Re + (1 - a)||V^||| 2 < (|H|^ 2 + \\v 2 f H .)\\vfv + T^zPW%- 

at L — a 1 — a 

Taking the L 2 inner product of (|3.2[) 9 with f , we are led to 
|-|(We||,||| 2 )T 



(3.4) 



\\ 2 = 2a(D(v)\f)-We((v-W)n+g a (f,Wv 1 )+g a (T 2 ,Vv)\f) 
S Aa 2 

< jprnb + —\\nh + CWe (HiJlUeUVnlUa + !|VC|| L2 ||r 2 j| L ^)||r|| L2 
+CWe ||Vui||ioo||f||| 2 

< i\\Vv\\b + ^f-\\nh + CWe HVtJH^dlnllfl, + \\r 2 \\ H ,)\\f\\ L 2 
+C*We ||Vui||^||t||£ 2 

< l\\Vv\\h + + l|Tl!l « 2 + l|r2|l ^ 2) + H V »iM We U*> 



for some d > and all S > 0. It follows from (O and (O that, for all t G [0, T], 



< 



Rc||S||| 2 + We||f||| 2 +/ ((l-a)||Vu||£ a + 2||r||£ a )<k 



TC 2 (1 + ||«i||£» (ifa) + \\v2\\U {m) )(\\v\\U(L*) + + s i m»ds 



+dj (^(l + \\r4m + \\r 2 \\ 2 H2 ) + 



\\y Vl \\ H 2)(We \\?\\l 2 )ds 



< TC 2 (l + 2B 1 )(\\v\\ 2 L ^ {L2) + \\e\\ 2 L ^ L2} ) + 5 I \\Vv\\ 2 L2 ds 



+Ci I (^(l + ||n|||r2 + ||r2||^) + ||V«i||H=)(We ||r||£ 2 )d S , 
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for some C2 > 0. Gronwall's inequality implies then 

Re ||«||i« (£a) + We ||f ||^ (i2) + f ((l-a)\\Vu\\ 2 L2 +2\\f\\ 2 L2 )d s 







< [rC 2 (l + 2B 1 )(||«||ioo (ia) + ||fl||io. (i3) ) + * / \\yv\\hds] 

Jo 

x[l + Ci(y(l + Hnlli.^) + ||r 2 ||^ 2) ) + VfWVvtW^^) 
x exp (Ci( T (l + Hnlli.^) + |N|^ (i?2) ) + v^UVwill^^)))] 



< [TC 2 (1 + 2B 1 )(|H| 2 ^ (L2) + ||0||i„ (L2) ) + S f \\Vv\\hd*\ 





x [1 + d(j (1 + 2B 2 ) + v^ 7 ^) exp (Ci(y(l + 25?) + y/TBl))], 

Setting (5 = min{Re , We , l-a}(4+8Ci cxp(2Ci))- 1 and T* = min{Ti, j-^gj, ^, 4 c 2 (1+2^ i+2§I cxp^d )) } ■ 
we see that for all T < T*, 

+ + J (||Vu||| 2 + ||f||| 2 ) ds < ^(\\vf LOO(L2) + lieilic^,) + J ||V«||i a ds). 

Hence, <E> is contractive as a mapping from Y(T*) to y (T+). The assertion of Proposition ^, ll thus follows 
from Lemma 13.21 

□ 

4. Proof of the Main Theorem 

Let (u,p,r) be the local solution to equation 11.11 constructed in Proposition 13. II We recall from this 
proposition that 

u G L 2 (0,T + ;# 3 ) n C{[Q,T i ];D{A)) with u' G L 2 (0, V) n C([0, T*]; L 2 ) and 
r G C([0,T*];iT 2 ) with r' G C([0, T*]; ff 1 ). 

Our proof for the existence of a unique, global solution to (|1.1|) is based on the following a priori 
estimates for u,t,u' and t'. 

Let us begin with an a priori estimate for r. To this end, we take the inner product of (|l.l[h with r 
and obtain 



(4.1) __|| T ||2 2 + || T ||2 2 = 2a (D(u)\T)-We (g a (r, V«)|r) < 2a||Vu|| i 2j|T|| i 2+CWe |!V M || ff2 ||r||| 2 . 



We d 

Similarly. 
We d 



||Vt||£ 2 + ||Vr|| 2 2 = 2q(V£>(u)|Vt) - We (V 5a (r, Vw)| Vr) - We {dm* <9 fe T y \d lT l - 



2 dr 

< 2a||V 2 u|| L 2||Vr|| i 2 + CWe ||Vu|| H 2||Vr|| 2 2 + CWe || V 2 m|| L 2 ||t|| l » || Vt|| L 2, 
and for i, j, k, I, m G {1, 2, 3} 

^^||V 2 r|| 2 2 + ||V 2 t||| 2 = 2a(V 2 D(u)\V 2 T) - We (d lm g a (T, Vu) <J '|8k ro T«) - We (d /mU fc d fc T^' |ft m r«) 

-We (a m u fc a fci T^|a im r^) - We (9^ fe 9 fcm r y |9 im r«) 
< 2a||V 3 w|| L 2||V 2 r|| i 2 + CWe ||Vw|| H 2||V 2 r|| 2 2 + 

CWe ||V 3 u|| ia ||T|| L =||V 2 r|| L2 + CWe || V 2 u\\ m || Vr|| ffl || V 2 r|| L 2 . 
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Combing the above three inequalities, we obtain 

\r\\ 2 H 2 + \\r\\ 2 H2 < 2a||Vu||fl a ||r|| fla +CWe||V«|| H »||r||| ra 



We d „ 



2 dt 

I a 



l\\r\\l 2 +Ca 2 \\^u\\% 2 +C^\\r\\% 2 , 



d " '12 , IUII2 ^ ^,2||V7„.I|2 ■ ~ We 2 



and thus 

(4.2) We ^||r||^ + ||r||^ < Cd'\\Vuf H , + C^\\r\\% 2 . 
Before estimating u, let us first apply the Helmholz projection P to (|l.l[h . This yields 

(4.3) Re (d t u + P((u • V)u)) + (1 - a)Au = Pdivr. 

Next, we estimate the term || Vm||#2 appearing in the right hand side of (|4.2[) . Corollary |2.2l and Remark 
12.31 imply that 

(4.4) !|V 2 U || ffl <C{\\Au\\ H , + \\Vu\\ L 2). 
We deduce from equation f)4.3[) that 

||VA u || i2 < J£_||Vut|U» + 1 ^-||VP((u-V)«)|| La + T ^— ||VPdivr|| ia 
1 — a 1— a 1 — a 

(4.5) < -5f_||Vut|U» + ^-||(«-V)«|| irI + -^— ||VPdivr|| £a 
By the Gagliardo-Nirenberg as well as by Sobolev's inequality, we have 

(4.6) iMii- < q|«||i 8 ||v«||t, < CHVuH^llVull^, 

which allows us to bound the term || (u • V)u||#i as 

\\(u-v)u\\ H , = ||v((«-v)«)|U a + ||(«-v)«|U a 

< ||Vu|||« + ||M|| L -i|V 2 u|| L 2 + ||u|| L 6||Vu|| L 3 

< C (\\Vuf H1 + |Mli,||Vu|| I 6 ||V 2 u || i2 + HVullillVull^) 

< C {\\Vu\\ 2 m + IIVull^HVull J, + HVullLllVullIx) 

(4.7) < C\\Vu\\ 2 H1 . 
Recalling Corollary 12.21 and Remark 12.31 we infer that 

(4.8) ||Vu|| H i <C(\\Au\\ L2 + \\Vu\\ L2 ). 
Combing now the estimates (|4.4 jl -(|4.7 jl with (|4.8|l yields 

||V U || 2 ff2 < C(p U || 2 2 + ||V U || 2 2 + ^^||V Ut || 2 2 + ^^||VPdivr|| 2 2 

^ +^\\M\% + J ^\\Vu\U 

Finally, estimate (|4.2p combined with estimate (|4.9p implies that 

We |||r||| 2 + Hrll^ + 2\\Vuf H2 < Kl {\\Au\\%, + ||V U || 2 2 + ||V U( || 2 2 + || VPdivr|| 2 2 ) 
(4.10) +C(\\Au\\i 2 + \\r\\% 2 + \\Vu\\t 2 ), 

for some Ki > 0. 
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Next, taking the inner product of (|4.3|) with u, we obtain 

^-^NIi» + (l-«)l|V«||i a = (divr|u). 

Adding this equation to equation (|4.ip . integrating by parts and using the fact that r is symmetric, 
yields 

^(Re||u|| 2 L2 + ^||r|| 2 L2 ) + (l- a )||V W ||i 2 + ^||r||i 2 = ~(j9 a (r, V«)|r) 

CWe „ „„ . 
< -2^-||T||L<»||V«|| L 2||r|| L2 

l-a„_ „, C*We 2 



which means that 



(4.H) |(He ||«||£ a + ^l|r|| 2 L2 ) + (1 - a)||V U || 2 2 + i||r|| 2 2 < ( f W ^ 2 IM| 4 g2 ■ 

Taking the inner product of (|4.3p with Am yields 

^IIV* + (1 - a)\\Auf L2 < ^\\Auf L , + -^||Pdivr|| 2 L2 + ^-||P(u • V)u\\\ 2 . 
A at A 1 — a 1 — a 

Further, since 

(4.12) ||P(u.V)u||£ a <C||(u-V)u||£ a < CNl£»||Vu||i. <C\\Vu\\l4Vu\\ m , 
it follows that 

(4.13) ^(Re ||Vu||i 2 ) + (1 - a)\\Au\\ 2 L2 < -^-||Pdivr|| 2 L2 + e||Vu||^ + P** ^ \\Vu\\ 6 L2 . 
at 1 — a e[i — ay 

Similarly, taking the inner product of (|4.3|) with u', and using (|4.12|) once more, we see that 



(4.14) -((1 - «)||V«||i 2 ) + Re ||u'||£ 2 < _|jPdivr|j£ 2 + £ ||Vu||^ + — ^— ||Vu||^ 
for e > 0. In view of (|4.8[> . (|4.13p and (|4.14p and by choosing e small enough, we are led to 



d_ 

dt 

(4.15) < K2 (||V U || 2 2 + ||PdivT|| 2 2 )+C||Vu||i 2 , 



((2Re + 1 - a)||Vuj| 2 2 ) + Re \\d t u\\ 2 L 2 + (1 - a)\\Au\\ 2 L 2 



for some k 2 > 0. Next, differentiating equations (|1 . 1[) i and with respect to t, and taking the inner 

product of the resulting equations with dtu and <9 t T, respectively, we obtain 

^~^\\ d M\h + (! - «)l|V«t||i a = (divrtlfttt) - Re ((«( • V)«|ft«), 



as well as 

We d 



\\dtT\\ 2 L2 + \\d t T\\ 2 L2 = 2a(D(u t )\d t T) - We ((u t ■ V)r|9 t r) 



2 dt 1 

-We (ff (r t) V«)|r t ) - We (<? (t, Vu t )\d t r). 



It follows that 
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1 d ( ,.„..„ We 



Re \\d t u\\i, + -7^\\dtr\\h ) + (1 - a)||V«t||l a + ^ll^lfo 



We n We 

< C\\Vu\\ H 2(Re \\d t u\\ 2 L2 + —\\d t r\\ 2 L2 ) + — \\d t u\\ L e\\Vr\\ L s\\d t T\\ L 2 

Za Za 

+ C^\\Vu t \\ L 4T\\ L ~\\d t T\\ L2 

za 

We o We 

< C\\Vu\\ H ,(Re ||a tU || 2 L2 + — ||a t r||i 2 ) + C— ||V Ut || L2 ||r|| ff2 ||5 t r|| L3 , 

za za 



and by Young's inequality that 



^ Re \\ dtU f L2 + ™\\ dtr f L2 ) + (i _ aJHVntlli, + -j|d t r|| 2 2 



dt V 11 2a 



Q 



r* Wp 2 fw P 2 

(4.16) < ellVulllp + -(Re 2 \\d t u\\h + ^\\d t r\\U + dMMft*. 

e 4or a [1 — a) 

Next, following an idea of Molinet and Talhouk [TB], we estimate Pdivr and curldivr, which will be 
then used in order to control ||Pdivr|| jji . In order to do so, we take the divergence of (jl.ip ^ and, using 
the incompressible condition, we obtain 

We (divTj + div ((u ■ V)r)) + divr + We div g a (r, Vu) = aAu, 

which, together with the equation of u, implies that 
(4.17) 

1 — a 1 1 — a 

We divT t + -divr = We [div ((u ■ V)r) + div g a (r, Vu) + Re [d t u + (u ■ V)it] + Vp. 

a a a 

Applying the Helmholtz projection P to this equation yields 

— -We PdivT t + -Pdivr = Re <9 t u + P|Re (u • V)u - ^— ^We (div((u • V)r) + div g a (r, Vu)) 1 . 
a a a 

Taking the inner product of the above equation with Pdiv r and integrating by parts, we deduce that 

l_a We_ d || Pdivr ||2 + — ||PdivT||? 2 = -Re (Vu t |r)+Re (F(u • V)w|Pdiv r) 
a 2 dt a 

1 — a 1 — a 

(4.18) We (Pdiv((u • V)r)|Pdivr) We (Pdiv 5o (r, Vw)| Pdivr). 

a a 

Note that 

Re |(P((w V)u)|Pdivr)| < Re ||uj| i6 || Vit|| L 3j|r|| ff i < CRe ||Vu|| L2 ||Vw|| H i||r|| H i 

£„_ „o CRe 2 



< ^\\Vu\\ 2 H1 + —- llVulli.llrll^, 



and 
1 — a 1 



1 




a 




a 




1 




a 




a 




1 




a 


a 



< We ||V U |M|r|| 2 H2 < £ -\\Wu\\ 2 H1 + — (1 , " } ||r|| 4 g2 , 

a o ear 

where we have used the Gagliardo-Nirenberg inequality (|4.6[) again. 

Moreover, a direct calculation yields 

—We |(Pdiv. 9a (r, V U )|Pdivr)| < ^||V«||^i + — 2 ~ ^ ||rj| 4 g2 . 
a o ear 
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which implies 
(4.19) 

^(i_^ W e ||Pdivr||| a ) + ^||Pdivr|| 2 2 < Re 2 ||V Ut ||| 2 + ||r||| 2 + e\\Vu\\ 2 H1 + C £ (||V U || 4 L2 + ||r||^) . 

Applying the curl operator to the equation (|4.17p we obtain 

1 - a . 1 

We curl div r t H curldivr = 

a a 

I - a 

Re [curl u t + curl ((« • V)w)] We curl div ((u ■ V)r) + curl div g a (T, Vit) . 

a 

Taking the inner product of this equation with curl div r yields 

1 — a We d „ , ,. .., 1 „ , ,. ,,, 

^TT + curldivr £ 2 + - curldivr £ 2 < 

a 2 at a 

-^-11 curl div r II | 2 + Ca ^ e \\Vu t \\ 2 L 2 + Re (curl((u • V)u)\ curldivr) - 
2a 2 

1 — a 1 — Gt 
(4.20) We (curl div ((u ■ V)r)|curldivr) We (curl div g a (r, Vit)|curl div r). 

a a 

Noting that 

[curl div ((u • V)r)]y = (u • V)(curldivr)y + (d u u k d k T ij - d jlU k d k T il ) 

+ (d iU k d kl T ij - d lU k d kj T U ) + ■<)!<,■■<),, 7'* - d ]U k d M T a ), 

we obtain 

— -We |(curldiv((u- V)r)|curldivr)| < C— -We (||V 2 d| L6 j|Vr|| L3 + II Vd| L ~ || V 2 rj| L 2) ||V 2 rj| L 2 
a a 

e..„ „, C(l - a) 2 We 2 „ n4 

EOT 

Moreover, 

Re |(curl((wV)u) | curldivr) | < CRe (|| |Vw| 2 || L 2 + || |w||V 2 w||| L 2)|| V 2 r|| L2 

< CRe (HVulU-HVullia + ||u|U6||V 2 u|U3)||V 2 r|U 2 

< CRe ||Vu||^||Vu|U=||r|U 2 < |||V U || 2 H2 + ^-HVull^lHI^, 

and a direct computation yields 

—We [(curl div <? a (r, Vu)|curl div r)| < E -\\Vuf m + C(1 ~ - ||r||^ 2 . 

a o ea z 

Substituting these estimates in (|4.20[) . we obtain 

(4.21) 

—We |||curldivr|| 2 2 + l||curldivr|| 2 2 < CaRe 2 ||Vu t || 2 2 + £ ||V U j| 2 ff2 + C 6 (||Vu||i, + ||r|| 4 ff2 ). 
a at a 

Putting together and (jOTjl yields 

4 ( l—^We (||Pdivr|| 2 2 + ||curldivr|| 2 2 )^ + - (||Pdivr|| 2 2 + j|curldivr|| 2 2 ) 
at \ a J a 

(4.22) < 2£||Vu|| 2 f2 + K3 ||Vu f |li2 + IM|i2+C e (||V U || 

for some K3 > 0. On the other hand, multiplying equation (|4.15p with "j^ 1 and adding to (|4.10p yields 

^(Wellrll 2 ,. + (K1 + ^ +1 - a) llWHi.) + ^±^||^||1 2 + + ||rj| 2 g2 + 2\\Vu\\% 2 

at 1 — a 1 — a 

(4.23) < k 4 (||V«||| 2 + UVutllia + l|Pdivr|| 2 ffl ) + C (\\Au\\ 4 Ll + ||r||^ + ||V«||i 2 + ||V U || 6 L2 ) , 
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for some K4 > 0. 

Finally, we estimate ||Pdivr||jji in the right hand side of (|4.23[) . Notice first that in view of the 
Helmholtz decomposition, we verify that curldivr = curlPdivr. Moreover, since div(Pdivr) = and 
(Pdivr) ■ v = 0, in virtue of Proposition 12. 51 there exists a constant Co such that 

(4.24) ||PdivT||^i < Co (||Pdivr||^2 + || curl div r|| ^ 2 ) . 

Then multiplying (|4.22p with a(«4Co + 1) and adding to (|4.23|) implies that 

d ^We Hi-Hip + (Kl + 1)(2Re +1 ~ a \ \Vu\\ 2 L2 + (1 - a)( K4 C + l)We (||Pdivr||i 2 + ||curldivr||2 2 ; 



dt V 1 - a 

(«i + !) Re ,| 9tU ||2 2 + p u ||2 2 + || r ||^ 2 + 2 j|v U ||^ 2 + ||PdivT||| a + ||curldivr||| 2 



1 - a 

< n 5 (e\\Vu\\ 2 H2 + \\Vu\\ 2 L2 + \\Vu t \\ 2 L2 + \\t\\ 2 l2 ) + C z {\\Au\\\ 2 + \\t\\ 4 h2 + \\Vu\\t 2 + ||Vu||£ 2 ) , 
for some K5 > 0. 

The term ||Vu t ||^ 2 in the right hand side of above can be absorbed into the left hand side by means 
of (|4.16p . Indeed, multiplying (|4.16p by yj^ 1 , adding the resulting equation to (|4.25p and choosing 

l-q 

k 5 (2-q) + 1 ' 



s = . /n °\ M , we infer that 



I (We |M|^ + (Kl + 1)( i me Q +1 - a) ||V,||| 2 + ^±^\\d t u\\h + ^^ll^lll, + 

(1 - a)( K4 C + l)We (||Pdivr||| 2 + ||curldivr||| 2 )) + Re ( Kl + ^ ||Q tU ||| 2 + 

1 — a 

\\Au\\ 2 L2 + \\r\\ 2 H2 + \\Vu\\ 2 H2 + HVutHi, + ; ! £f ± 4ll^T|l!= + l|Pdivr||| 2 + ||curldivr|| 2 L2 
< C {\\Au\\t 2 + ||r|| 4 ff2 + ||V«||i a + HVttllk + \\d t u\\ 4 L2 + ||a t r||l 2 ) + n e {\\Vu\\ 2 l2 + \\t\\ 2 l2 ) , 

for some kq > 0. Finally, using the assumption < a < 1 and multiplying (|4.1ip with £s±I anc l adding 
it to (f4725|) yields 

d ( B£^±1) || U ||2 2 + ^^||r||», + We | M |i, 2 + (Kl + 1)( i me Q +1 - a) ||V,||| 2 

+ Re (K 5 + l) ^ + We («» + 1) ]]dtTf a){K4Co + 1)We ( || Pdivr || 2 + 

1 — a 2a{L — a) 

HcurldivrlliO) + Re(Kl + 1) ||^ll| 2 + \\M\h + \\r\\h + W^Wh + W^Wb 
1 — a 

+ { H +1 \ \\dtT\\h + ||Vu|| 2 L2 + ||r||| 2 + ||Pdivr||| 2 + !|curldivr||| 2 
a(l — a) 

(4.25) < C {\\Au\\i 2 + \\t\\% 2 + ||V U j| 4 2 + ||V U || 6 L2 + ||ft«||* a + ||a t r||| 2 ) . 

After all these calculations, we are glad to define the functions F : [0, 00) — >• R, G : [0, 00) — > R and 
H : [0, 00) -> R as 

F(t) := (1 - a)(K 4 C + l)We (||Pdivr||| 2 + ||curldivr||^ 2 ) 
Re (« 6 + l) 2 We (« 6 + l). 



\l 2 



2a{l - a) 



li 2 +We ||r||^ 



- (/tl + 1)( i 2Re Q +1 - Q) l|y M ||| 2 + Re 1 (K5 + 1) ||^||j 2 
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and 

;= Re ( Kl + 1) 2 + 2 + 2 + a + a 
1 — a 

+ i^±ll| M | 2 + || V „||3 a + ||r||| 2 + ||Pdivr||l 2 + ||curldivr||| 2 , 
a(l — a) 

and 

H(t) := ||9 t n||| 2 + \\Au\\h + \\rf m + \\d t r\\ 2 L2 + ||Vu||£ a + l|Vu||£ 2 . 
We now rewrite (|4.25p as an inequality of the form 

(4.26) j t F{t) + G{t) < CH(t)G(t) 

and estimate H(t) in terms of F(t). To this end, observe that by (|4.3p . (|4.12p and (|4.8p we deduce that 

(4.27) \\Au\\ 2 L2 < C (\\Vu\\l 2 + \\d t u\\ 2 L2 + ||Pdivr|| 2 L2 + ||V U || 6 L2 ) . 
Hence, there exists a constant Mi = A/i(Re , We , a) > such that 

(4.28) H(t) < Mi(F{t) + F(t) 2 + F(tf), t > 0. 
Substituting (|438]) into ([426]) . we get 

(4.29) ^ J F(t) + (l-CA/i(F(t)+F(t) 2 + f(t) 3 ))G(t)<0, * > 0. 

We are now finally in the position to estimates F(i). In order to do so, define S > small enough, such 
that Sq + S 2 + Sq < 2cTKh' wnere C i s the constant appearing in (|4.29[) . 

Assume that the differential inequality (|4.29p holds for all t > and F being absolutely continuous 
and G being nonnegative. Then 

(4.30) F(t) < 5 for all t > provided F(0) < 5 - 

Assume that this assertion were not true. Let t\ > be the first time where F(t) > So. Then 

(4.31) F(h)=So and F(t) < S for all 0<t<h. 
Consequently, for all < t < t\, 

I - CMi(F(t) + F(t) 2 + F(tf) > 1 - CAhiSo + S 2 + S 3 ) > ^ 
Assertion (|4.29[) implies now that 

(4.32) ^i F ^ + \ Gi ^ ~ for a11 
Integrating (|4.32p from to t\. we obtain 

(4.33) F(h) + I f 1 G{s)ds < F(0) < S , 

* Jo 

which contradicts (|QT|l . Thus (|OP)l holds true. 

Given this fact, the proof of Theorem II. 1 1 can now be finished easily. In fact, let T* be the lifespan of 
the local solution (u,p,r) given in Proposition 13. 11 Assuming that (|1.5p holds with Eo to be determined 
below, we verify that 

F(0) < C(\\u Q \\ 2 Hl + \\t \\ 2 h2 + K(0)||| 2 + ||r t (0)||| a ) < C(e 2 + e A ). 

Choose now £o such that C{e\ + Eq) < So- Then, in virtue of (|4.30[) and the proof of (|4.33[) . we obtain 

(4.34) sup F(t) + i / G(t)dt < S . 
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In particular, it follows from the definition of F(t) and G(t), ([4~34]) . (|4~27ll and (|4~8|) that 
sup (||u(t)||i, (A) + \\u'(t)\\ 2 L2 + \\r(t)\\ 2 H2 + !|r'(t)|| 2 L2 ) 

<t<T* \ y ' / 



0<t<T* 

(llVu^H^ + \\Vu\t)\\l 2 + \\ T (t)\\h + \\r(t)\\h)dt < C. 

We thus deduce that the local solution (u,p,r) can be extended for all positive times. This completes 
the proof of Theorem 11.11 □ 
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